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Abstract 

We present a method for the propagation of partially coherent radiation using coherent mode 
decomposition and wavefront propagation. The radiation field is decomposed into a sum of inde- 
pendent coherent modes. Each mode is then propagated separately using conventional wavefront 
propagation techniques. The summation of these modes in the plane of observation gives the co- 
herence properties of the propagated radiation. As an example, we analyze propagation of partially 
coherent radiation transmitted through a single circular aperture. 

PACS numbers: 42.25.Kb, 42.25.Bs, 41.60.Cr 



1 



INTRODUCTION 



With the construction of the third generation synchrotron sources (see for example pQ) 
partially coherent radiation in the x-ray range has become available. The advent of the 
free-electron lasers [21 E] has tremendously increased the amount of coherent flux, that can 
be used by experimentalists. Techniques like coherent x-ray diffractive imaging (CXDI) [U- 
[6] explicitly utilize the high coherent flux and promise new insights in structural biology, 
condensed matter physics, magnetism and other correlated systems [7]. 

Careful planing of the coherence-based experiments is required, which means that the 
knowledge of the beam properties at the experimental station is in high demand. Although 
many computational methods have been developed to calculate the beam profile at the 
sample position, most of them do not provide the coherence properties of the beam in the 
plane of observation. In addition, the majority of these calculations are based on the ray 
tracing, or the wave propagation approach and are valid only in the limit of incoherent or 
fully coherent radiation, respectively. At the same time, partially coherent radiation that is 
mostly available nowadays at the 3rd generation synchrotron sources and free electron lasers 
is not covered by these methods. 

Here, we present a general approach, which can be applied to partially coherent wavefields, 
and which is capable of calculating both the intensity profile of the beam as well as the 
transverse coherence properties of the radiation at any position in the beamline. It is based 
on the results of statistical optics, where the radiation field is described in terms of correlation 
functions. In this work we represent the correlation function of the radiation field as a sum of 
coherent modes [8] . Each mode is propagated separately through the beamline optics using 
a wave propagation technique. At the sample position the correlation function is given as a 
sum of the propagated modes. The intensity profile and the transverse coherence properties 
of the radiation can be readily extracted from this correlation function. A similar approach 
was used in the frame of geometrical optics [H] and for advanced phase retrieval methods in 
coherent imaging |10j . 

We applied this technique to calculate the radiation properties of a Gaussian Schell-model 
source. As an example, source parameters of the free electron laser in Hamburg (FLASH) 
were considered and propagation of the partially coherent radiation through a pinhole was 
analyzed. 
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METHODS 



Definitions 

In the theory of partial coherence the mutual coherence function, r(ri, r 2 ; T), plays the 
dominant role. It describes the correlation between two complex scalar values of the electric 
field, E(r; t), at different points ri, r 2 and at different times t and t + r [8j [TT] 

r(r 1 ,r 2 ;r) = ( J E(r 1 ;t)*E(r 2 ;t + r)) T , (1) 

where the brackets (• • • )t denote the time average. It is assumed that the averaging is 
performed over times T that are much longer than the fluctuation time of the field |T2]. 
When we consider the propagation of the correlation functions in space, it is convenient to 
introduce the cross spectral density function, W(rx, r 2 ; u), as the Fourier transform of the 
mutual coherence function [8 J 

W( ri , r 2 ; u) = J r(n, r 2 ; r)e^dr, (2) 

where u is the angular frequency of the radiation. By definition, when the two points ri and 
r 2 coincide, the cross spectral density represents the spectral density [13] of the radiation 
field, I(t;lo) = W(r,r;u). The normalized cross spectral density is known as the spectral 
degree of coherence 

W{r u r 2 ;u) 

= . ==, < \fi 12 \ < 1. (3) 

^/(rijwji^ju) 

The modulus of the spectral degree of coherence is often measured in interference experi- 
ments as the contrast of the interference fringes. To characterize the transverse coherence 
properties of the wavefield by a single number the normalized degree of coherence can be 
introduced [2] 

= /iWW^fdr^ (4) 
(J/(r;c)dr) 2 

Mode decomposition and propagation of the correlation function 

It can be shown that under very general conditions the cross spectral density can be 
decomposed into a sum of statistically independent coherent modes [S] 

W(r lt T 2 ;u) = Y,Pn{u)E* n {r x -u)E n {r 2 -u), (5) 



where j3 n (u) are the eigenvalues and E n (r; uS) are the eigenf unctions of the integral equation 



In particular, the mode decomposition ^ can be applied to planar secondary sources [8], 
where the cross spectral density [15], W(ri,r 2 ) = W(ui,U2',Zo), of the radiation field is 
given in the source plane with the transverse coordinates u = (x,y). The coordinate z is 
defined along the optical axis and the position of the source is at z = 0. 

The mode decomposition is convenient in the analysis of the propagation of partially 
coherent radiation, when only a small number of modes is required to describe the cross 
spectral density. The coherent modes can be propagated separately along the optical axis 
and the cross spectral density function, W(ui,u 2 ; z), can be calculated at any position z 
along the optical axis via eq. ^ replacing the modes E n (u,Zo) by the propagated modes 
E n (u, z). The eigenvalues, (3 n , remain unchanged during the propagation. 

The propagation of the field for each mode over a distance z along the optical axis can 
be performed by a wave propagation technique. In the case of free space propagation this 
can be done utilizing the Huygens-Fresnel principle [16] 



where E n (u', zq) is the wavefield of the mode in the source plane and E n (u, z) is the propa- 
gated mode at position z. The propagator P z (u, u') is given by 



where k = u/c is the wavevector, r is the distance between the points (u',z ) and (u, z), 9 
is the angle between the line joining (u', zq) to (u, z) and the optical axis, and x(9) is the 
obliquity factor with x(0) = 1 and < x{@) ^ 1- 111 the paraxial approximation, when small 
angles 9 are considered, the propagator reduces to 



In general, the propagation of the partially coherent radiation through an arbitrary ar- 
rangement of the optical elements in a beamline can be performed in the following steps 

1. Decomposition of the cross spectral density, VF(ui, u 2 ; z ), of the source into coherent 






(7) 





modes E n (u, z ) according to eq. 
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FIG. 1: Partially coherent radiation is generated at the source. The circular aperture is positioned 
25 m downstream from the source. The coherence properties of the radiation are analyzed 45 m 
downstream from the aperture. 

2. Propagation of all modes E n (u, zq) from the source plane through the optical elements 
to the observation plane using a wave propagation technique. For example, equation 
([7]) can be used to propagate all modes from the source to the first optical element 
at the position z. For thin optical elements the transmitted modes are given by 
E out (u, z) = T(u)E(u, z), where T(u) is the transmission function. At the next step 
the transmitted modes are propagated to the next optical element using eq. (|7]). 

3. Finally, after performing the previous step for all optical elements present in the beam- 
line each mode is calculated in the plane of observation, and the cross spectral density, 
W(ui, U2; z), is determined by eq. 



Gaussian Schell-model 

A useful model to describe the radiation properties of partially coherent sources is the 
Gaussian Schell-model [8]. This model has been applied for the analysis of the radiation 
field generated by optical lasers [TT] , third generation synchrotron sources [H] and x-ray 
free-electron lasers (T8H2Q]. In this model the cross spectral density in the source plane, 
W(ui,u 2 ; z ), is given by 

W(u 1: u 2 ; z ) = exp — — — exp — — , (8) 



v 2al 2a* ) *\ 2g 2e y , 

where a x , a y is the source size and £ x , £ y is the transverse coherence length of the source in 
the horizontal (x) and vertical (y) direction, respectively. Due to the symmetry of the Gaus- 
sian Schell-model the total cross spectral density at the source factorizes into the horizontal 
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FIG. 2: (top) The total intensity distribution at the source position, (bottom) The intensity 
distribution of the nine lowest modes at the source position. The length of the white scale bar is 
150 fim. 

and vertical components 

W(ui,u 2 ;zo) = W(x 1 ,x 2 ;z ) ■ W(y 1 ,y 2 ; z ). (9) 

The modes E n and their corresponding contributions /3 n can be found for each direction 
separately 

W(xi, x 2 ; z ) = ^ PnK( x u z )E n (x 2 ; z ), 

n 

and a similar expression is valid for the y direction. The analytical solution of the integral 
equation ^ for the Gaussian Schell-model in each direction is known in the form of the 
Gaussian Hermite-modes. The eigenvalues in this model are described by a power law [S]. 
The total cross spectral density is given by 

VK(ui,u 2 ; z ) = (u 1 ;z )E nm (u 2 ;z ), (10) 

n,m 
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FIG. 3: (top) The total intensity distribution behind a circular aperture with (a) 3mm and (b) 1 
mm diameter, (bottom) The intensity of the lowest nine modes behind the aperture. The length 
of the scale bar is 2 mm. 

where ^^(u; z ) = E n (x; zq) ■ E m (y; zq) are the well known laser resonator TEM nm modes 
[2 X j and (3 nm = PnflL. The total radiated intensity can be calculated by 



I{u;z ) = W(u,u;z ) = ^ f3 nm \E nm (u; z )\ 2 . 



(11) 



We want to note here that our approach is not limited to the Gaussian Schell-model. If the 
cross spectral density of the field, W(u 1 ,u 2 ; z ), at the source is known, the corresponding 
modes E n and their contributions (3 n can be calculated from the integral equation ^ (see 
for example Ref. [22])- 



PROPAGATION OF PARTIALLY COHERENT RADIATION THROUGH A CIR- 
CULAR APERTURE 

We simulated the propagation of the partially coherent radiation generated by a Gaussian 
Schell-model source through a circular aperture (see Figure [I]). As source parameters we 
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FIG. 4: (top) The total intensity distribution in the plane of observation calculated for a circular 
aperture, (a) 3mm and (b) 1 mm in diameter, (bottom) The intensity distribution of the lowest 
nine modes in the plane of observation. The length of the scale bar is 5 mm. 

used values reported for FLASH |2] operating at a wavelength of 13.5 nm: the full width 
at half maximum (FWHM) source size is 160/zm and the FWHM divergence of the beam is 
90 /irad. In the frame of the Gaussian Schell-model (see eq. ([8])) these values correspond 
to a x = o y = 68 /im and £ x = £ y — 62 /mi |18j . For these parameters the contribution of 
each mode at the source position is given by the power law f3 nm = 0.41 n+m . Modes with 
a contribution of less than 1 % were neglected in our simulations. In total, 21 modes with 
n + m < 5, where n = 0, 1, ... 5, m — 0, 1, ... 5, were used in the calculations presented 
here. The intensity distribution of the total field and the nine lowest modes at the source 
position are shown in Figure [2} 

We considered the following geometry in our simulations: the pinhole is positioned 25 m 
downstream from the source and the radiation properties are analyzed 70 m downstream 
from the source (see Figure [TJ. Such an arrangement is typical for experiments performed in 
the unfocused beam at FLASH. We applied the general scheme of propagation of partially 
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coherent radiation described earlier for this experimental geometry. The propagation in the 
free space was performed using equation ^ with the propagator in the paraxial approxima- 
tion. Different pinholes with the diameter from 5 mm to 1 mm were analyzed. The pinhole 
transmission function, T(u), was defined as 



T(u) 



1 for |u| < d/2 
elsewhere 



where d is the diameter of the pinhole. It was convolved with a 200 fim wide (FWHM) 
Gaussian function to account for imperfections of the pinhole edges. The intensity distribu- 
tion of the total beam and the lowest modes behind 3 mm and 1 mm pinholes are presented 
in Figure |3j Figure [4] shows the same in the plane of observation. 

One readily sees in Figure [3] (a) that in the case of the 3 mm pinhole the first modes, which 
dominate the radiation field, are just slightly affected by the aperture. The 1 mm pinhole, 
however, substantially cuts all modes, including the dominant ones (see Figure [3] (b)). The 
intensity distribution of each mode in the plane of observation is similar to the intensity 
distribution leaving the pinhole. Additional intensity modulations due to the scattering on 
the edges of the aperture are observed in Figure [3] (a). This can be attributed to the Fresnel 
diffraction effects, which are stronger for sharper pinhole boundaries. The Fresnel number, 
d 2 1 (XL), where A is the wavelength of the radiation and L is the distance from the pinhole to 
the detector, is 15 for the 3 mm pinhole and 1.6 for the 1 mm pinhole in this experimental 
geometry. We want to note here that due to Fresnel diffraction, small variations in the 
propagation distances might change these intensity modulations significantly. 

Finally, the cross spectral density, W(ui, u 2 ; z), was determined in the plane of observa- 
tion. We present here results for the horizontal direction 

W(xi, x 2 ; z) = f3 nm E* nm (x l , y x = 0; z)E nm (x 2 , y 2 = 0; z). (12) 

nm 

Due to the symmetry of our scattering geometry, the same result is obtained in the vertical 
direction. The modulus of the cross spectral density, |W(a;i function of the 

transverse positions, x\ and x 2 , is shown in Figure [5] (a-d) for different apertures. The 
modulus of the spectral degree of coherence |/i(Ax)| = |//(— Ax/2, Ax/2) | as a function of 
the separation Ax for the same apertures is presented in Figure [5] (e-h) (red solid line). This 
calculation corresponds to the measurements of the contrast in a double pinhole experiment 
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FIG. 5: Transverse coherence properties of the radiation in the observation plane, (a-d) The 
modulus of the cross spectral density \W(xi,X2] z)\. (e-h) The modulus of the spectral degree of 
coherence |//(Ax)| as a function of the separation Ax. Simulations were performed with the pinhole 
diameters of (a,e) 5 mm, (b,f) 3 mm, (c,g) 2 mm, and (d,h) 1 mm. The insets in (e-h) show the 
intensity distribution I(x) in the horizontal direction. Red solid lines show the simulations with 
the presence of the circular aperture. The blue dashed lines show the same functions obtained 
without the pinhole. 

with varying pinhole separation Ax and the center of the double pinhole positioned at the 
optical axis of the beam. For comparison, the spectral degree of coherence, |/i(Ax)|, for 
the same geometry and source parameters, but without the aperture is shown by the blue 
dashed line. The intensity profiles, I(x), for the same aperture sizes are shown in the insets 
of Figure [5] (e-h) . 

As a result of our simulations, we notice that the 5 mm pinhole does not affect the 
transmitted radiation. The intensity profile I(x) as well as the modulus of the spectral 
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degree of coherence |/i(Aa;)| in the observation plane calculated with and without the pinhole 
are the same. For the smaller pinhole diameters of 3 mm and 2 mm the size of the beam 
decreases, but the modulus of the spectral degree of coherence is not significantly altered. 
Our simulations suggest, that in the present geometry down to the pinhole sizes of 2 mm 
no significant changes in the coherence properties of the beam in the observation plane are 
expected. Only for the smallest pinhole size of 1 mm the values of |/i(Ax)| are significantly 
enhanced at large separations (see Figure [5] (h)). In this case the spectral degree of coherence 
cannot be described by a single Gaussian function (compare to results of Ref. [23J). However, 
we should note here that these separations are much larger than the beam size and will be 
difficult to access in a real experiment. 

We analyzed as well the normalized degree of coherence £ ^ and the available photon 
flux P behind each pinhole (see Table 1). The transmitted photon flux was calculated by 
the equation P = f J oni (u)du/ f / in (u)du, where /j n (u) and I «i( u ) are the intensity (see 



eq. (11)) incident on and behind the aperture. If no pinhole is present in the beamline, 
then obviously the photon flux P = 100 % and the normalized degree of coherence has the 
value determined by the source parameters ( = 18 %. Results from the Table 1 show that 
the normalized degree of coherence, (, is significantly increased for the smaller pinholes. 
However, this happens at a loss of the transmitted photon flux P. It is interesting to note 
that the product P ■ £, which may be considered as the amount of the coherent photon flux, 
is about 20 % after transmission through the larger pinholes. It drops down to a value of 
only about 10 % for the 1 mm pinhole. In the latter case almost a fully coherent beam , 
C = 78 %, is achieved with 10 % of the transmitted radiation. 



TABLE I: The transmitted photon flux, P, and the normalized degree of coherence, behind the 
aperture. Four different pinhole diameters are analyzed. 



no pinhole 


P -- 


= 100% c = 


18% 


5 mm 


P 


= 97% C = 


19% 


3 mm 


P 


= 69% C = 


29% 


2 mm 


P 


= 39% C = 


46% 


1 mm 


P 


= 10% c = 


78% 
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CONCLUSIONS 



In conclusion, we have presented a computational method, that allows the calculation of 
the transverse coherence properties as well as the beam intensity profile of partially coherent 
radiation at any position in the beamline. Our approach can be easily implemented, since it is 
based on the wave front propagation, for which several powerful computational methods are 
already developed. The important extension to the conventional wave propagation methods 
is the consideration of all contributing modes. 

We have applied this method to describe the propagation of partially coherent radiation 
through a circular aperture. Parameters typical for FLASH were used in the simulations. 
Our analysis shows, that the presence of the pinhole does not increase the transverse co- 
herence length substantially, even if the beam is significantly cut by the aperture. At the 
same time, the normalized degree of coherence is significantly increased for smaller pinholes 
at the expense of the available photon flux. A straightforward extension to our approach 
would be its application to a beamline with a number of optical elements. Our simulations 
suggest that a careful calculation is important to have a realistic picture of the beam prop- 
erties in the observation plane for a concrete realization of the optics in the beamline. An 
interesting question, for example, would be how imperfections of the optical elements affect 
the transverse coherence length of the radiation. 

This approach can be effectively used to describe the radiation at free electron lasers, 
which is predominantly coherent and contains a small number of statistically independent 
modes. It can be also applied to third generation synchrotron sources, which are highly 
coherent in the vertical direction and less coherent in the horizontal direction |14j . However, 
in this case a significantly larger number of modes has to be taken into account. 
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